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Owing to its numerical simplicity, a two-dimensional two-electron model atom, with each electron 
moving in one direction, is an ideal system to study non-perturbatively a fully correlated atom 
exposed to a laser field. Frequently made assumptions, such as the "single active electron" -approach 
and calculational approximations, e.g. time dependent density functional theory or (semi-) classical 
techniques, can be tested. In this paper we examine the multiphoton short pulse-regime. We 
observe "non-sequential" ionization, i.e. double ionization at lower field strengths as expected from 
a sequential, single active electron-point of view. Since we find non-sequential ionization also in 
purely classical simulations, we are able to clarify the mechanism behind this effect in terms of 
single particle trajectories. 

PACS Number(s): 32.80.Rm 



I. INTRODUCTION 

Several theoretical approaches were able to reproduce experimentally observed ion yields in multi-electron ionization, 
at least qualitatively (see e.g. [^]). Most of them are based on a "single active electron" (SAE) point of view [^,||. 
A new impact on the research in this field had the discovery of the so-called "knee" or "shoulder" in the ionization 
yields of helium exposed to a laser pulse Q . This means that double ionization occurs many orders of magnitude 
more strongly at intensities where, according to a sequential SAE scenario, almost no He++ should be present. Early 
after the experimental observation of this non-sequential ionization (NSI), two possible mechanisms were suggested in 
order to explain it. Corkum proposed a rescattering scenario [|[ where the first electron revisits the core and ionizes 
the second electron collisionally. Fittinghoff et. al. suggested a "shake off" effect where the second electron ionizes 
due to the sudden loss of screening of the core by the first electron. Walker et. al. |^ concluded by analyzing their 
experimental data that a rescattering process is not able to explain the observed yields. Their arguments are based 
on the absence of a rigorous threshold in the He++ yields. Instead they propose "that NSI occurs via a simultaneous 
two-electron ejection either through a shake off or threshold mechanism involving some form of electron correlation" . 
Recently, the NSI mechanism has been clarified within the intense-field many-body S'-matrix theory |^ . It was shown 
"that the dominant mechanism behind the observed large probability of laser-induced double escape is a quantum 
mechanical process of absorption of photon energy by one of the electrons which is shared cooperatively with the 
other electron through the Coulomb correlation^^ This mechanism for the NSI process was independently deduced 
from ID He-studies where the model atom had been exposed to a low frequency, short pulse laser field ||]: "[...] 
before the outer electron disappears completely, the inner electron is already sufficiently strongly excited so that it 
leaves the atom within a short time interval later. It is during this time interval that the correlated double ionization 
takes place." Simulations where the outer electron is calculated in the SAE way but the inner one feels (in a second 
computer run) the time dependent potential created by the outer one, succeeded in reproducing the NSI- "knee" 
This result is also a strong indication that the suggested mechanisms as quoted above are, indeed, the correct ones. 
However, there is no detailed physical picture how this energy sharing between the outer and the inner electron takes 
place. 

Our calculations were performed for a relatively high frequency {uj = 0.4 a.u.) and a very short pulse duration (6 
optical cycles) while in ^ a low frequency short pulse was used. Since we are rather in the multiphoton-regime than 
in the tunneling domain, the occurrence of NSI might be surprising at all. Indeed, in our calculations NSI is relatively 
weak compared to the many orders of magnitude effect for ionization of helium in strong low frequency laser light. 
However, with the help of our additional classical simulations we are able to provide (i) a detailed physical picture 
how NSI takes place in terms of one-particle trajectories, and (ii) a proof that NSI, in its essence, is not a quantum 
mechanical effect. 
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Because the full quantum mechanical numerical simulation of helium exposed to a laser field is an extremely 
demanding task [|l^, approximate approaches are desirable. Among these, Hartree-Fock- [O-^, time-dependent 
density functional- (TDFT) and semi-classical molecular dynamics-calculations |20H23] are most frequently 

used. Especially the latter method succeeded in reproducing the "knee" [^. On one hand the molecular dynamics 
calculations are very appealing and instructive since particle trajectories and the single particle energies can be traced. 
On the other hand the additional "Heisenberg-force" which must be introduced in order to avoid instabilities where 
one electron falls into the "black hole" (i.e. the nucleus) while the other one ionizes, is somewhat artificial and may 
evoke objections against the results produced by this method. 

The time-dependent Hartree-Fock method was found to be problematic in the framework of multi photon ionization 
]T3|-p!5[. Results from TDFT, in principal an exact approach, depend on the choice of the effective exchange-correlation 
potential [Q. Another disadvantage of this procedure is that only the total electron density n(r,t) = \(pi{r,t)\'^ 
is calculated and the single particle orbitals (fiir^t) are physically meaningless in a rigorous sense. 

The study of systems where the motion of each electron is reduced to one spatial dimension has a relatively long 
tradition. Potentials of the form —Zj\Jx^ + e, so-called "soft core Coulomb potentials", provide an energetic Rydberg- 
like scaling |2^] and lead to results, qualitatively similar to those from full 3D calculations. Two ID electrons are a 
two-dimensional system which is tractable with nowadays computers. Two ID electron-systems are used to study non- 
perturbatively autoionization [|6|, ionization of a negative ion validity of time-dependent Hartree-Fock theory 
for the multiphoton ionization of atoms [|l3| -p!5|, and, most recently, two-electron effects in harmonic generation and 
ionization [||. 

This paper is organized as follows. In Section || the model system is introduced. In Section HI results from the 2D 
quantum calculations are presented. Section |^ is devoted to an analysis of the results in terms of a SAE-approach. 
In Section |^ we present the results from a time-dependent density functional-calculation and in Section VI we discuss 
our classic al p article simulations within which the NSI scenario can be clarified. Finally, we summarize and conclude 
in Section Ivil. 



II. THE ID HELIUM MODEL 



The two ID electrons with coordinates x and y interact with the core and with each other through a "soft core"- 
interaction, i.e. + e and 1/ ^{x — yY + e, respectively, and with the field E{t) through the dipole term 

{x + y)E{t) (atomic units (a.u.) will be used throughout this paper). Thus the total Hamiltonian reads 

.r/ N 1 1 2 2 1 . x^. X 

H{x,y,t) = --—--—- + + X + y)E(t). 1 

2dx^ 2dy^ V^^+e vV + e ^/{x^^y)^Te 

The desired ground state energy can be tuned by varying e. We used 

e = 0.55 

in our calculations which leads to the ground state energy 

£o = -2.897 a.u. 

on our numerical grid. Sq is approximately the ground state energy for the real 3D helium atom which is —2.902. 
One may prefer thinking in terms of one 2D particle which moves in the somewhat peculiar 2D potential 

2 2 1 

V{x, y,t) = — -== -== + = + (x + y)E{t) (2) 

yx'^^e V^^ + e ^J{x-yY^e 

instead of the two electrons interacting with each other. The potential (H) and the ground state energy level are 
shown in Figure |l| for the three constant fields E = 0.0, 0.1, and 0.616. The electric field E tilts the fieldfree potential 
around the axis y = —x. 

In order to estimate at what fieldstrengths strong single and double ionization should occur, it is advantageous to 
calculate the classical critical fields. However, the commonly used method of equating the initial ground state energy 
level to the maximum of the barrier formed by the atomic potential and the external field leads to an unphysically 
small critical field ^^^^ = 0.009 for single ionization. In our numerical simulations in the next sections we will find a 
negligible ionization probability for such low field strengths. Like in the case of hydrogen-like ions, the electron is not 
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able to move beyond the point where the energy level touches the lowest point of the energetic barrier |2^] . Therefore 
the "real" critical field is higher. On the other hand, total ionization should be strong at latest when the ground state 
energy level exceeds the electron-electron repulsion ridge in x = ^/-direction. This is the case at = 0.616 (see 

Figure ||c). 

III. QUANTUM CALCULATION 

We used a spectral method to solve the time-dependent Schrodinger equation (TDSE) 

d 

with the Hamiltonian (||). In order to keep the numerical effort as small as possible we chose a rather high laser 
frequency 



and a very short pulse covering 6 cycles, 



uo = 0.4 



T = 6— = 94.248 . 

CO 



The pulse envelope was sin -shaped, thus 

E{t) = E sin^ {^t) sin ujt (3) 

for < t < T. With the frequency chosen and E not greater than 1 we are in the multiphoton-domain since the 
Keldysh-parameter y^\£o\/{2Up) is not much less than unity over the whole intensity region of interest. Up is the 
ponderomotive potential £^^/(4c<;^), i.e. the mean quiver energy of a free electron in the laser field. A spatial grid 
spacing Ax = Ay = 0.4 and a temporal one At = 0.1 was found to be sufficient. We chose the grid size large enough 
so that almost no probability density reaches the boundaries within the pulse time T. The ground state wave function 
was determined by propagating a Gaussian seed function in imaginary time. A contour plot of the ground state 
probability density is shown in Figure ^ Its energy is Sq = —2.897. 

If we assume that one electron is already ionized we are left with the ID version of hydrogen-like helium IIe+. The 
ground state energy of this system, 

t) = f --^ - ^= ) t), (4) 



dt ' ' ' \ 2dx^ / 



^x- 



was determined to be f + = —1.920. Therefore, the ionization energy for the "outer" electron is Eq — = —0.977. 
The removal energy for the outer electron in real 3D He is 0.9 a.u. 

The probability density y, t)p during the pulse for the peak field strength E = 0.3 is shown in Figure ||. One 
looks perpendicularly from above onto the illuminated, logarithmically scaled probability density. 

Density fiowed along the ±x- or ±?/-axis corresponds to single ionization while double ionization happens when 
both \x\ and \y\ are significantly greater than the width of the ground state. We use this simple picture to define 
our ionization probabilities at the end of the laser pulse for total, double and single ionization P, and P+, 

respectively. 



1 



/a pa 
/ \'^{x,y,T)\^dxdy, (5) 
-a J —a 



/oo na 
/ \^{x,y,T)\^dxdy 
-ooJ —a 

/a nOO 
/ \<i{x,y,T)\''dxdy 
-aJ —OO 

^{x,y,T)\^dxdy, (6) 
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(7) 



We chose a = 5. In Figure ^ the integration areas corresponding to P, P++ and P+ are shown. 

Of particular interest is the evolution of the probability density when the double ionization regions \x\ > 5 and 
\y\ > 5 are occupied for the first time. In Figures ^ and || these time intervals are shown for E = 0.2 and 0.7, 
respectively. The delay between the snapshots is a quarter optical cycle in each case. 

Let us first analyze the E = 0.2-case in Figure ^. Till n = 2.75 optical cycles mainly single ionization has occured, 
i.e. the probability density is still located along the axes. At time n = 3.00 two density jets enter the region x^y > 
(indicated by two arrows). They are clearly separated by the Coulomb repulsion ridge along y = x. These density 
jets represent states where both electrons are on the same side of the nucleus. Therefore the Coulomb repulsion is 
relatively high and the jets tend to flow back towards the single ionization channels. This reflux of probability density 
is supported by the electric field which has its maximum a quarter of a cycle later, at n = 3.25. However, some density 
passes the single ionization channels and appears in the regions x > O^y < and x < 0, > 0, respectively (see arrows 
in the n = 3.25- and in the n = 3.5-plot). Now, the two electrons represented by this density are on opposite sides of 
the nucleus, and both are ionized. So we conclude that although the electric field amplitude E is not sufficiently strong 
to double ionize by simply tilting the xy-plane. Coulomb repulsion and electric field can lead to double ionization by 
acting together constructively. Therefore we propose Coulomb repulsion assisted laser acceleration to be responsible 
for NSI. 

At higher intensities {E = 0.7, Figure ^) we also observe the two density jets entering regions of the x?/-plane where 
X and y have equal signs (see n = 2.00 till n — 2.25). However, reflux of these jets is obviously not essential to 
stimulate subsequent emission of probability density into regions x > 0, < and x < 0, > 0, respectively. This 
can be seen in plot n = 2.25 where already bursts of density leave the single ionization channels although the jets are 
not flown back yet. Here, double ionization is mainly due to the strongly tilted x?/-plane. 

Note that these backflowing jets which support double ionization are of course absent in any SAE model while they 
are expected to be included in TDDFT since there both electrons are allowed to respond to the field simultaneously. 
Furthermore, purely classical simulations should show a similar NSI scenario as the quantum density current in the 
xy-plane does, i.e., electrons which collide while moving in the same direction (the two jets) followed by a subsequent 
turn of one electron which crosses the origin and finally leaves in the opposite direction. All these presumptions will 
be confirmed in the next sections. 

We performed several runs varying the peak fieldstrength E. In Figure ^ the ionization degrees P, P+, and P++ 
are plotted vs. the intensity I = E^. There are also shown the ionization probabilities Pgequ and P • Pgequ which 
results from solving (^) in the external laser field. As long as P++ is small P • Pgequ would be the probability for He++ 
production if double ionization occurs purely sequential, i.e., the first electron does not interact with the residual ion 
anymore and the second electron remains non-excited after the emission of the first one. Obviously this completely 
sequential scenario cannot explain the correct ionization degree P++. Double ionization occurs earlier than predicted 
by P • Pgequ- This corresponds to the experimentally observed "knee" in the He ion yields, although our He++-curve 
does not show such a pronounced shape due to the fact that we are in the multiphoton and not in the tunneling 
regime as in We have also performed some runs using a; = 0.2 and saw an increasing deviation from the SAE 
He++ curves, i.e. an increasing "knee" . The Pgequ and P • Pgequ curves cross the correct P++ result since depletion of 
the He+ ions is not taken into account. 



IV. SINGLE ACTIVE ELECTRON ANALYSIS 

In the previous section we have already compared the correct ionization degree P++ with the product P • Pgequ for 
the rate when the second electron ionizes from its He+ ground state without interacting with the electron already 
freed. The ionization degree P • Pgequ was found to be too small in the intensity region of interest 0.02 < I < 0.07, as 
depicted in Figure 0. 

Now, we may try to describe the ionization probability P+ in the single active electron approach. In order to do 
this we have to find an appropriate Zeff and to solve the equation 

The first electron has an ionization potential So — = —0.977. We found that Zeff = 1.117 yields such a binding 
energy. In Figure H the resulting ion yields P+ are shown and compared with the exact result P+. The ionization 
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degree calculated from (||) is too high over the whole region where mainly single ionization takes place. As soon as 
double ionization occurs the curve differs also qualitatively from the exact one: there is a dip in the exact P and P+ 
curves around / = 0.03 which is absent in the P+ result. We conclude that the second electron shares some energy 
with the escaping first electron which leads to a decreased single ionization probability. Since in a similar work for 
a very low frequency [|j the SAE curves are found to fit well, the deviations we observe might be mainly due to the 
relatively high frequency we have chosen. For higher frequency laser fields the "cracking up" of the initial ground 
state wave function into an "inner" and an "outer" electron may be too slow to be well described by an SAE ansatz 
for the "outer" electron. 

One may object that although tuning Zeff in order to fit the binding energy of the "outer" electron leads to an 
over-estimation of the ion yields, there might be a certain combination of Zeff and an effective soft-Coulomb-parameter 
eeff which does both: providing the correct binding energy and the right ionization probabilities. However, we tried 
Zeff = 1 and eeff = 0.398 which leads also to the desired binding energy —0.977. The resulting curve for the single 
ionization yields are also shown in Figure ||. It overestimates the single ionization too. 

V. DENSITY FUNCTIONAL THEORY 

Density functional theory (DFT) is a powerful tool in determing multielectron atomic structures (see e.g. [|l6|,0 for 
an overview). It has been shown that a Hohenberg-Kohn-type theorem exists also for time-dependent phenomena. 
Therefore the existence of an effective potential which transforms the problem of N interacting electrons to that of 
N non-interacting ones is proved. The non-interacting electrons move in an effective potential which is a functional 
of the total electron density only. The problem reduces to finding an appropriate effective potential which includes 
all relevant exchange and correlation effects. 

There is some recent work on the field of time-dependent density functional theory (TDDFT) applied to laser 
ionization of atoms 

In the case of a singlet two-electron-system like our model atom of helium, there is only one occupied Kohn-Sham- 
orbital (p{x^t). The total electron density is 

n{x,t) = 2\(p{x,t)\^. (9) 
There are no exchange contributions, and neglecting correlation effects leads to the time-dependent Hartree equation 

X — X' r + e 



Hartree's independent electron-model was used by Geltman already in 1985 to analyze experimental results in the 
multiple ionization of atoms js^Js^. 

The ground state energy we obtained by solving ([lo| ) in imaginary time is = —2.878. 

In Figure ^the comparison between the TDDFT-results and the exact ones are made. The ionization yields for He+ 
and He++ are observables, of course. Thus they are functionals of the density n(x, t) and, due to the simple relation 
(H), explicit functionals of the Kohn-Sham-orbital density \^{x^t)\^. We want do adopt the simple "integration-over- 
areas" picture in order to calculate the ionization, as described in Section |l|. According to Ullrich et. al. [|l^ we 
proceed, for the time being, as follows: with 



Pks 



:= 1 - / \Lp{x,t)\^dx (11) 

J —a 



the probability for neutral helium is the product of the probabilities for each orbital to be non-ionized. Thus 

P°s = (1 - Pks)^ (12) 

For the single and double ionization 

P+s = 2Pks(1 - Pks), (13) 
P^^ = Pis (14) 
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follows and the total ionization clearly is 

^KS=l-^KS- (15) 

The three curves corresponding to Pks ' ^ks shown in Figure ^. The agreement is quite bad. The total 

and single ionization degree are overestimated as in the SAE calculation. 

However, the total ionization Pks fits well the exact P if one avoids assigning any physical relevance to the Kohn- 
Sham-orbital t) and proceeds instead as follows: Pks as defined in (^l|) is the probability to find any of the two 
electrons outside the interval [—a, a] since the physical total electron density is n{x^t) = 2\(f{x^t)\'^. Therefore the 
probability for total ionization should be simply 

^^KS = ^'ks, (16) 

and 

^KS = 1 - ^KS. 

^KS according (|l6|) is also depicted in Figure ^ The agreement between Pks and P is excellent. The dip around 
/ = 0.03 is well reproduced. Since the dip is absent in the SAE results, the onset of NSI seems to be included in the 
TDDFT although only the simple Hartree effective potential is taken. 

There is no simple way to deduce Pj^g and Pks^ without claiming physical significance of the Kohn- Sham-orbital 
(f{x^t) although the existence of purely density dependent functionals Pj|g [n] and Pk^ [n] are proved by the Hohenberg- 
Kohn-theorem. Eqs. (p!2|-|l5|) would be valid if the correct wave function ^(x, t) was the product of the Kohn-Sham- 
orbitals, (p{x^t)(p{y^t). A plot of \(p{x^T)\'^\(p{y^T)\'^ for E = 0.3 is shown in Figure |lo|. Clearly, there is a gridlike 
pattern imprinted due to the construction of the wavefunction as a pure product. This leads to a totally different 
angular distribution in the x?/-plane. A similar behaviour is observed in time-dependent unrestricted Hartree-Fock 
calculations [|l^. 

However, for higher field strengths both electrons are ionized rapidly and subsequently behave as free and almost 
independent electrons. Thus the total wave function should develop a grid-like pattern. This is shown in Figure pT| 
where the peak field strength E = 0.7 was chosen. The TDDFT result is also plotted. 



VI. CLASSICAL SIMULATIONS 

We solved the classical equations of motion according the Hamiltonian (|^) and the electric field (||) for a micro- 
canonical ensemble of the two electrons. We traced the one-particle-energies 

12 1 

^x{x, y) = -x^ , + =, (17) 

^y{^. y) = If - + ' (18) 

^ V r + e V(x -yy 



The total energy is 



S{x,y,t) = £x{x,y) -^£y{x,y) = + {x^y)E{t). 

^/{x-yy +e 



Each electron is considered to be ionized when Sx{x^y) > or £y{x^y) > 0, respectively. There is no unique way in 
defining the single particle energies (|l^-|l^ since the electron-electron term may be shared between the two electrons 
in various ways |2^. However, this has little influence on the ionization degrees since in the case of single or double 
ionization the distance between the two electrons is normally large at the end of the pulse. 

The initial conditions were chosen to meet the quantum mechanical ground state energy So = —2.897. Fortunately, 
the resulting ion yields were not sensitive to the choice of the ensemble of initial conditions. Instead of taking several 
initial positions and momenta we started with one "mother" -configuration at t = and vary the time ton where the 
laser pulse sets in. We tried several mother-configurations to ensure insensibility of the resulting ion yields. 

We would like to mention that a classical treatment of an ID model helium has also been undertaken in in the 
framework of stabilization of multielectron atoms. 
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The results are shown in Figure |T2|. The single ionization is strongly overestimated. This is due to the fact that the 
ionizing electron gains energy at the expense of the still bound electron which occupies a state of quantum mechanically 
forbidden low energy, i.e. an energy below —1.920. This is shown in Figure 13 where the two one-particle-energies £x 
and £y are plotted for a representative single ionization event at £^ = 0.1. 

This behavior could be prevented if one introduces a velocity dependent "Heisenberg" -potential However, 
for our purpose to study the NSI-mechanism this is not necessary. We have calculated also the classical SAE single 
ionization process in the potential V{x) = —Zq^I\Jx^ + e with Z^^ = 1.117 as in Section IV. The resulting yields 
are lower and show a rapid increase at the critical intensity which is 0.05 for the potential used. This is the typical 
behavior in pure classical simulations. Note that double ionization is already observed where according the SAE 
approach even no classical single ionization occurs. 

The result of the SAE calculation He+ ^He++ was multiplied with the probability for He+-production from the 
full 2D run. Below / = 0.25, the classical critical intensity for sequential He++-production the probability vanishes, 
as expected. Thus the intensity region 0.03 < I < 0.2 is the classical NSI regime we are particularly interested in. 
We examined each double ionization event in that region. The ionization dynamics of two representative examples in 
terms of single-particle-energies and trajectories is depicted in Figures |lj and |l5|. 

For the purpose of comparing our classical results with the quantum mechanical probability density current in 
the xy-pldbue we look at the electron trajectories in Figure |l5|. One observes that the electron which leaves first is 
accompanied by the other electron moving in the same direction (corresponding to density flowing into regions of the 
xy-pldine where x and y have equal sign). The Coulomb interaction is strong within that half cycle. Then, as the 
first electron ionizes, the second one turns around and moves in the opposite direction (corresponding to probability 
density passing one of the x?/-plane's axes). Thus, the second electron in NSI leaves the atom approximately half 
laser cycle after the first one. The dynamics of the classical test particles corresponds to the temporal and spatial 
evolution of the quantum probability density as described in Section III. 

In Figure [l^ we show a representative example for double ionization at a higher field strength {E = 1.0). At this 
field strength the sequential pathway is more probable than NSI. The temporal delay between the ionization of the 
two electrons is greater (1-3 cycles) but still small since our pulse is ramped strongly over 3 cycles only. In longer 
pulses the temporal delay between the ejection of the two electrons would be even greater. 



VII. CONCLUSIONS 



In this paper we have confirmed the recently proposed mechanism for the NSI process, namely the ionization of the 
second electron by Coulomb-interacting with the outer partner. We have traced the NSI scenario by observing the 
quantum mechanical probability density as it evolves in space and time. We have identified the Coulomb repulsion 
assisted laser acceleration of the inner electron to be responsible for NSI. Our classical simulations have supported this 
point of view and contributed a detailed picture of how NSI happens in terms of one particle energies and trajectories. 
Moreover, we have shown that "non-sequential" means "within half an optical cycle" and that NSI is not an essentially 
quantum mechanical effect. We showed that in the frequency and pulse duration regime under consideration, SAE 
ionization yields are in poor agreement with the exact results. TDDFT reproduces the total ionization probability 
very well, even for the simplest effective potential available, namely the Hartree-potential. TDDFT fails if one claims 
physical relevance for the Kohn- Sham-orbit als by separately constructing single and double ionization yields. However, 
since generally only the knowledge of the total electron density is necessary, e.g. for the calculation of high harmonics 
generation, TDDFT should produce very accurate results [|l7[ . 



ACKNOWLEDGEMENT 



The author would like to thank P. Mulser and R. Schneider for helpful discussions. This work has been supported 
by the European Commission through the TMR Network SILASI (Super Intense Laser Pulse-Solid Interaction), 
No. ERBFMRX-CT96-0043. 



7 



* http:/ /www. physik.th-darmstadt.de/tqe/ 
[1] S. Augst, D. D. Meyerhofer, D. Strickland, and S. L. Chin, J. Opt Soc. Am. B 8 858 (1991) 
[2] K. J. Schafer, Baorui Yang, L. F. DiMauro, and K. C. Kulander, Phys. Rev. Lett. 70 1599 (1993) 

[3] Baorui Yang, K. J. Schafer, B. Walker, K. C. Kulander, P. Agostini, and L. F. DiMauro, Phys. Rev. Lett. 71 3770 (1993) 
[4] D. N. Fittinghoff, P. R. Bolton, B. Chang, and K. C. Kulander, Phys. Rev. Lett. 69 2642 (1992) 
[5] P. B. Corkum, Phys. Rev. Lett. 71 1994 (1993) 

[6] B. Walker, B. Sheehy, L. F. DiMauro, P. Agostini, K. J. Schafer, and K. C. Kulander, Phys. Rev. Lett. 73 1227 (1994) 
[7] Andreas Becker and Farhad H. M. Faisal, J. Phys. B: At. Mol. Opt. Phys. 29 L197 (1996) 

[8] D. G. Lappas, A. Sanpera, J. B. Watson, K. Burnett, P. L. Knight, R. Grobe and J. H. Eberly, J. Phys. B: At. Mol. Opt. 
Phys. 29 L619 (1996) 

[9] J. B. Watson, A. Sanpera, D. G. Lappas, P. L. Knight and K. Burnett in: Multiphoton Processes 1996, Inst. Phys. Conf. 
Ser. No 154, ed. by P. Lambropoulos and H. Walther, pp. 132, lOP Publishing Bristol 1997 
[10] J. Parker, K. T. Taylor, C. W. Clark and S. Blodgett-Ford, J. Phys. B: At. Mol. Opt. Phys. 29 L33 (1996) 
[11] K. C. Kulander, Phys. Rev. A 36 2726 (1987) 
[12] K. C. Kulander, Phys. Rev. A 38 778 (1988) 

[13] M. S. Pindzola, D. C. Griffin, C. Bottcher, Phys. Rev. Lett. 66 2305 (1991) 
[14] M. S. Pindzola, P. Gavras, and T. W. Gorczyca, Phys. Rev. A 51 3999 (1995) 
[15] M. S. Pindzola, F. Robicheaux, and P. Gavras, Phys. Rev. A 55 1307 (1997) 

[16] E. K. U. Gross, J. F. Dobson, and M. Petersilka, Density Functional Theory of Time- Dependent Phenomena in: Topics in 

Current Chemistry Yo\. 181 pp. 81, Springer Berlin, Heidelberg, 1996 
[17] S. Erhard and E. K. U. Gross in: Multiphoton Processes 1996, Inst. Phys. Conf. Ser. No 154, ed. by P. Lambropoulos and 

H. Walther, pp. 37, lOP Publishing Bristol 1997 
[18] C. A. Ullrich, S. Erhard and E. K. U. Gross, Density- Functional Approach to Atoms in Strong Laser Pulses in: Super 

Intense Laser Atom Physics IV, ed. by H. G. Muller, NATO ASI series 3/13, pp. 267, Kluwer 1996 
[19] C. A. Ullrich and E. K. U. Gross, to appear in Comments on Atomic and Molecular Physics 
[20] D. A. Wasson and S. E. Koonin, Phys. Rev. A 39 5676 (1989) 
[21] P. B. Lerner, K. J. LaGattuta, and J. S. Cohen, Phys. Rev. A, 49 R12 (1994) 
[22] P. B. Lerner, K. J. LaGattuta, and J. S. Cohen, J. Opt. Soc. Am. B, 13 96 (1996) 
[23] P. B. Lerner, K. J. LaGattuta, and J. S. Cohen, Laser Physics, 3 331 (1993) 
[24] C. A. Ullrich, U. J. Gossmann, and E. K. U. Gross, Phys. Rev. Lett. 74 872 (1995) 
[25] J. Javanainen, J. H. Eberly and Qichang Su, Phys. Rev. A 38 3430 (1988) 

[26] D. R. Schultz, C. Bottcher, D. H. Madison, J. L. Peacher, G. Buffington, M. S. Pindzola, T. W. Gorczyca, P. Gavras, D. 

C. Griffin, Phys. Rev. A 50 1348 (1994) 
[27] R. Grobe and J. H. Eberly, Phys. Rev. A 48 4664 (1993) 
[28] D. Bauer, Phys. Rev. A 55 2180 (1997) 

[29] M. D. Feit, J. A. Fleck, Jr., and A. Steiger, J. Comp. Phys. 47 412 (1982) 

[30] R. M. Dreizler, E. K. U. Gross, Density Functional Theory, Springer Berlin, 1990 

[31] Erich Runge and E. K. U. Gross, Phys. Rev. Lett. 52 997 (1984) 

[32] Sydney Geltman, Phys. Rev. Lett. 54 1909 (1985) 

[33] S. Geltman, J. Phys. B: At. Mol. Opt. Phys. 21 47 (1988) 

[34] Maciej Lewenstein, Kazimierz Rzazewski, and Pascal Salieres in: Super-intense laser-atom physics, ed. by B. Piraux, A. 
L'Huillier and K. Rzazewski, NATO ASI series B: physics Vol. 316, pp. 425, Plenum 1993 



FIG. 1. 

The 2D potential V{x,y) (|) for E = 0.0 (a), E = 0.1 (b) and E = 0.616 (c). At E = 0.1 the initial ground state 
level cuts the effective potential within the single ionization channels. For E = 0.616 the ground state level even 
exceeds the potential ridge along x = y. At latest at that field strength strong double ionization is expected. 

FIG. 2. 

Contour plot of the ground state probability density. The energy is So = —2.897. The electron-electron repulsion 
along X = y clearly leaves its fingerprint on the wave function (butterfly-shape). Such an asymmetric shape is absent 
in corresponding Hartree-Fock jlj] or simple DFT ground states. 



FIG. 3. 
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The probability density after n optical cycles for ^ = 0.3. One looks perpendicularly from above onto the illuminated 
xy-plane and the logarithmically scaled probability density 7/, Till n = 2 mainly single ionization takes place 
(probability density along the axes). Afterwards also regions l^/l > 5 are occupied which corresponds to double 
ionization. The Coulomb-repulsion ridge along y = x can be clearly identified at later times. 



FIG. 4. 

The areas to be integrated over in order to calculate the probabilities for (from left to right) total, double and single 
ionization are shaded. The parameter a was chosen to be 5 a.u. 

FIG. 5. 

The probability density for £^ = 0.2 at n = 2.5, 2.75, 3, 3.25, 3.5, 3.75 optical cycles. At time n = 3 two probability 
density jets are emitted into the region x,y > 0. Atn = 3.25 they are partly flown back and density appears in 
regions x < O^y > and y < O^x > 0. 

FIG. 6. 

The probability density for E = 0.7 at n = 1.5, 1.75, 2, 2.25, 2.5 optical cycles. Double ionization is mainly due to 
the strongly tilted X7/-plane. 

FIG. 7. 

The 2D calculation results for total, single and double ionization (bold, +). The SAE result for the He+ ^He++ 
is also drawn (O). Multiplication with P (see text) leads to the expected sequential double ionization degree (*) as 
long as P++ <^ P+. The Pgequ and P • Pgequ curves intersect the correct P++ result since depletion of the He+ ions 
is not taken into account. 



FIG. 8. 

SAE results for the outer electron compared with the exact 2D yields (bold, +). The curve plotted with connected 
* was calculated using Zgff = 1.117 and egff = e = 0.55. For the O-curve Zgff = 1, egff = 0.398 was chosen. In both 
cases the SAE ionization degrees are too high. 



FIG. 9. 

Comparison of the TDDFT results with the exact 2D solutions. The total ionization degree Pks (A) matches 
nicely the exact curve (bold, +). If one claims physical relevance of the Kohn-Sham orbitals one gets the O-curves 
for P^g* (total), Pj^g (single) and P^^ (double) ionization which agree poorly with the exact probabilities (see text 
for a discussion). 



FIG. 10. 

The TDDFT probability density at the end of the pulse for ^ = 0.3. The gridlike pattern is due to the construction 
as a pure product of the Kohn-Sham orbitals. There is no displacement of probability density along y = x, oi course. 



FIG. 11. 

The probability densities at the end of the laser pulse for E = 0.7. In plot (a) the exact density, in (b) the TDDFT 
result is presented. The exact density (a) develops a gridlike pattern since both electrons were ionized rapidly and 
subsequently behaved like almost free electrons, i.e., the wave function becomes more and more a pure product of 
single particle wave functions. However, the agreement with the product of the Kohn- Sham-orbit als (b) is poor even 
at those high fieldstrengths. 

FIG. 12. 
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The classical yields (bold, *) for Pd (total), (single) and P^~^ (double) ionization. The single ionization is 
strongly overestimated (compare with the exact quantum mechanical results, drawn dashed) due to the classical 
effect discussed in the text. The classical SAE results for the outer (Pjj", □) and the inner electron {Pj^~^ ^ A) are also 
plotted. There is classical NSI at / = 0.04 even when no sequential single ionization should occur! 

FIG. 13. 

A representative example for classical single ionization in terms of the single particle energies (|l^,|l8|). The inner 
electron (solid) drops below the quantum mechanical He+ binding energy —1.920. This leads to an overestimation 
of the single ionization probability since the outer electron (dashed) can gain more energy during Coulomb collisions 
(see peaks in the energy-curves) with the inner electron as allowed quantum mechanically. 

FIG. 14. 

Two representative classical NSI scenarios in terms of the single particle energies (|l7| , [l8|) . 

FIG. 15. 

The particle trajectories corresponding to Figure 0. The electrons become free within approximately | an optical 
cycle. Before both electrons ionize they move together in the same direction for a quarter of a cycle until one electron 
turns and finally vanishes in the opposite direction. This has to be compared with the quantum dynamics in Figure 

i 

FIG. 16. 

A representative example of classical sequential double ionization at E = 1.0. The temporal delay between the 
ejection of the two electrons is 3 half cycles and would be even greater in a more adiabatically ramped pulse. 
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Figure 1 ("A two-dimensional, two-electron ..." by D. Bauer) 
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Figure 2 ("A two-dimensional, two-electron ..." by D. Bauer) 
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Figure 4 ("A two-dimensional, two-electron ..." by D. Bauer) 
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Figure 5 ("A two-dimensional, two-electron ..." by 



D. Bauer) 
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Figure 6 ("A two-dimensional, two-electron ..." by D. Bauer) 
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Figure 7 ("A two-dimensional, two-electron ..." by D. Bauer) 




Figure 8 ("A two-dimensional, two-electron ..." by D. Bauer) 
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Figure 9 ("A two-dimensional, two-electron ..." by D. Bauer) 
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Figure 10 ("A two-dimensional, two-electron ..." by D. Bauer) 
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Figure 11 ("A two-dimensional, two-electron ..." by D. Bauer) 
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Figure 12 ("A two-dimensional, two-electron ..." by D. Bauer) 
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Figure 13 ("A two-dimensional, two-electron ..." by D. Bauer) 
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Figure 14 ("A two-dimensional, two-electron ..." by D. Bauer) 
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Figure 15 ("A two-dimensional, two-electron ..." by D. Bauer) 
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Figure 16 ("A two-dimensional, two-electron ..." by D. Bauer) 
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